Role of structural fluctuations in the insertion into complex host matrices 
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A coupling of structural and thermodynamic fluctuations in the course of various-type 
insertion processes is investigated within a combination of Gibbsian statistics and the 
information theory approach. It is shown that the coupling makes it possible to restore 
(at least partially) the information, inaccessible from experimental tests. This enables 
one to make physically reasonable predictions under limited information on the system. 



I. INTRODUCTION 



Intercalation processes find their application in 
many technologically important domains, such 
as the design of hydrogen-storage systems [1], 
rechargeable high-energy batteries, electrochromic 
devices, (see Ref. [2] for a review), electroactive 
polymers, and superconductors [3]. Microscopi- 
cally the insertion is a complex process involv- 
ing many effects (e.g. the charge transfer, elastic 
response of the matrix, the permselectivity, etc). 
Even well-characterized materials, like crystalline 
or layered compounds, exhibit rather complicated 
elastic properties, involving restructuring, staging 
and random distortion of the galleries [4]. For 
amorphous (or porous) matrices [5] the situation 
is complicated by the host heterogeneity, that may 
change in the course of intercalation. For instance, 
structural heterogeneity of electroactive polymers 
changes because of the polymer swelling. Complex 
geometry and poorly characterized energetics [6] 
make it practically impossible to construct a me- 
chanical model. Therefore, statistical mechanical 
investigations in this domain face difficulties, in- 
duced by the system complexity. In this situation 
it seems reasonable to develop a new theoretical 
scheme, that would be capable of making reason- 
able predictions under restricted information. 

For this purpose we combine the standard sta- 
tistical mechanics and the information theory ap- 
proach [7]. The latter has been introduced as an 
alternative (with respect to the standard statis- 
tical mechanics) tool for the statistical descrip- 
tion of many-body systems. Such an approach has 
been successfully applied to a description of non- 
equilibrium steady states and the self-organization 
[8] in various complex (physical, chemical, biolog- 
ical, etc.) systems. One of the main advantages 
of the information theory is its capability of de- 
scribing the equilibrium and nonequilibrium states 
in the framework of a unified scheme, resembling 
the usual Gibbs statistics, provided that a suit- 
able information (entropy) measure is chosen. On 
the other hand, the approach does not rely upon 
microscopic details. This allows one to study the 
systems deviating from the standard equilibrium 
conditions. Quenched systems, like fluids adsorbed 
in porous media [9], spin glasses [10,11] or hquid 



glassformers [12] could serve as relevant examples. 

In this paper we consider complex systems in 
which the information at a microscopic level is not 
sufficient for making a link between the system me- 
chanics (Hamiltonian) and its statistics (probabil- 
ity distribution) . Then the we can split the system 
into subsystems with different levels of description. 
One is a dynamic subsystem, which evolves accord- 
ing to a Hamiltonian, containing unknown param- 
eters. The latter reflect a coupling to a stochastic 
subsystem which is governed by a probability dis- 
tribution. This differs from the usual quenching 
in, at least, two aspects. The dynamic subsystem 
can influence the stochastic one and their coupling 
determines both the thermodynamic behavior and 
the shape of the probability distribution. For in- 
stance, by the analogy with the intercalation sys- 
tem [13-15], we can expect a change of the matrix 
volume upon the fluid absorption. Secondly, one 
necessarily deals with two entropic impacts: the 
thermodynamic entropy (due to the dynamic sub- 
systems) and the information entropy (due to the 
probability distribution) . Constructing from these 
two a suitable entropy measure, we investigate the 
strain distribution and the guest thermodynamics 
in homogeneous and heterogeneous matrices. 



II. MAXIMUM ENTROPY ESTIMATION 

Let us assume that the states of a system are la- 
belled by a continuous variable x, with p{x) being 
the probability of a state x. Nevertheless, because 
of the system complexity, the probabilities are not 
known a priori. The only available information is 
on a constrained quantity 



Q 



dxp{x)q{x) 



(1) 



which is the expectation value of a physical observ- 
able (e.g. energy, heat flux, etc.). In addition the 
probability distribution is assumed to be normal- 
ized 



dxp{x) = 1 



(2) 



Given this information can we estimate another 
relevant quantity G = J dxp{x)g{x)7 For this pur- 



pose one has to restore the probabihty distribution 

p{x). 

The basic role in this approach is played by the 
entropy S that gives a measure of missing informa- 
tion [7] concerning the system state. If wc have the 
entropy as a functional of the probability distribu- 
tion S = S\p{x)], then, according to the informa- 
tion theory approach [7,8], the probability can be 
estimated through the entropy maximization un- 
der constraints (1),(2). If we take the entropy (or 
the information) measure in the Shannon form 



S = — J dxp{x) lnp{x) 



(3) 



then such a variation procedure gives 

e-79(a:) 



p{x) 



Z = j da;e-T«(^) (4) 



where 7 is a Lagrange multiplier which has to be 

determined from the constraint (1). For instance, 
if Q is associated with the equilibrium internal en- 
ergy, then p{x) becomes the conventional Gibbs 
distribution, with 7 = l/{kT). 



III. COMPLEX SYSTEMS 

As is discussed above, by complexity we mean a 
situation when the available microscopic informa- 
tion on a system is not sufficient for its description 
in terms of a Hamiltonian. Thus, only one subsys- 
tem, say {s}, is governed by a Hamiltonian H[a, s], 
while the rest (the surrounding) is specified by a 
set of relevant parameters {cr} which appear with a 
probability distribution P{<t). Therefore, we have 
a coupling of a dynamic system to a stochastic one, 
which however influence each other. This means 
that P(cr) is unknown, the only available informa- 
tion concerns with an observable quantity 



j daP{a)E{a) 



(5) 



that gives, for instance, a characteristic energy 
scale for the surrounding. In the case of spin 
glasses the counterparts of {s} and {a} systems 
are the spin variables and the random fields (or 
exchange constants), respectively. For porous me- 
dia, {cr} can be identified with the matrix (e.g. 
pore size) and {s} - with the adsorbate degrees of 
freedom. 

For a given configuration {cr} one can calculate 
the partition function of the dynamic subsystem 



Z{a) 



j {ds)e 



-l3H[<7,s] 



(6) 



where /3 = l/{kT) is the inverse temperature. 

Then all the thermodynamic characteristics are 
known. For instance, the free energy 



F{a)=--\nZ{a) 



and the internal energy 



t/(a) = --lnZ(a) 



(7) 



(8) 



determine the thermodynamic entropy 

ST{cT)=f3[U{a)-F{a)] (9) 

Then the relevant thermodynamic quantities can 
be obtained by averaging over all realizations of a 
. For instance, the average free energy is given by 



F = (F(a)) = J {da)P{a)Fia 



) (10) 



Note that, in contrast to the conventional equilib- 
rium, now we have two entropies - the thermody- 
namic entropy 

St = j {da) P {a) St (a) (11) 

and the one related to the information on the prob- 
ability distribution 

Si = - j {da)P{(j) In P{a) (12) 

For a given distribution P{a) the value of St is 
the average thermodynamic entropy of the cou- 
pling {s}-{(t}. If, however, the distribution is not 
known, then St can be viewed as an additional en- 
tropy measure for inferring P{a) through the vari- 
ation procedure. 

In the case when the surrounding affects the dy- 
namic subsystem, but not vice versa (pure quench- 
ing) P((t) can be determined through the Sj maxi- 
mization under the constraint (5) and the normal- 
ization condition. This gives 



P(a) = 



-7B(a) 

Zo 



Zo = j dae-^^(") (13) 



where E{(j) can be estimated through model or 

even scaling arguments. For instance, E{cr) = aa^ 
is the energy to create a field of the magnitude 
cr. As is discussed above, the Lagrange multiplier 
7 should be determined from the condition (5). 
Once the probability distribution is determined, 
the thermodynamic quantities can be calculated 
straightforwardly. 

Such a decoupling of the two subsystems is how- 
ever an idealization, which could be more or less 
acceptable depending on a concrete situation. In 
general, the surrounding should respond to the 
evolution of the dynamic subsystem. For instance, 
the porous materials could change the volume 



upon accommodation of the fluid, like in the case 
of intercalation systems [13-15]. Also, it has been 
demonstrated [16] that the effect of the quench- 
ing on the thermodynamics of HCl-ice interfaces 
weakens with increasing ice film thickness. 

If the subsystems do influence each other, then 
the probability distribution must be evaluated 
starting from the total entropy 1^{St,Si), that 
describes the uncertainty concerning the state of 
the overall system. Obviously, E(0, Sj) = Sj and 
T,{St, 0) = St when there is no loss of information 
on one of the subsystems. The main problem is to 
construct Ti{SttSi) for nonvanishing St and Sj. 
For instance, we can formally expand around the 
perfectly ordered state 5/ = 0, 5t = 

E(6't, Si) = S(0, 0) + kSt + XSi + SStSi + ... 

(14) 

As a first approximation the total entropy 

'E{St,Si) can be estimated as a quasi-additivc 
combination (we drop I](0,0) as an irrelevant con- 
stant) 



'^{St, Si) = Si + kSt 



(15) 



where k is a parameter reflecting the coupling be- 
tween the dynamic and the stochastic subsystems. 
Then the distribution function can be inferred by 
maximizing (15) under the constraint (5). This 
leads to 



P{a) = 



Z = J dae-^'^^''^+''^-^ 



(16) 



IV. STRAIN DISTRIBUTION IN THE 
COURSE OF INSERTION 



In this section we analyze the interplay of these 
two impacts considering various aspects related to 
insertion systems. In this case the guest plays the 
role of s-subsystem, whose coupling to the host 
matrix is given by the matrix strain e which is a 
counterpart of cr-subsystcm. We arc interested in 
determining the strain distribution P{e) and the 
guest thermodynamics St, taking into account the 
host-guest coupling St{£)- 



A. Homogeneous matrix 



As an illustration of our approach we consider 
the simplest example of guest insertion into a " ho- 
mogeneous" host matrix. "Homogeneous" means 
that the matrix properties (the structure, the bind- 
ing site distribution, etc.) are similar for different 
matrix domains (e.g. graphite galleries). It is well- 
known that the matrix becomes strained with in- 
creasing guest concentration x. For a given x the 
actual strain is a sum of the domain strains e^. 
For homogeneous matrices we may expect that all 
the domains are equivalent, such that J2i =ne, 
where n is the number of domains. 

Therefore, the host-guest coupling leads to a 
strain distribution P{£) and affects the guest ther- 
modynamics. The strain-dependent guest entropy 
St{£) can be expanded around e = 



where the Lagrange multiplier 7' should be deter- 
mined from the constraint (5). It is known [17] 
that the entropy S't(c) determines the probability 
of thermodynamic fluctuations for a given value of 
a according to 



St{£) = St{0) + ane + ... 



(19) 



(17) 



It should be emphasized that P(p{cr) does not de- 
scribe the fluctuation of a, but the fluctuation of 
thermodynamic variables related to {s}-subsystcm 
(e.g. temperature, density, etc) for a given value 
of a. Then the resulting distribution (16) can be 
represented as 



Here St{0) is the guest entropy in the case when 
the matrix is not strained. This corresponds to the 
lattice gas model [2] describing the guest adsorp- 
tion on a rigid lattice. Note that all the quantities 
above {ST{0),a,e) are concentration dependent. 
The strain distribution P{e) can be determined by 
maximizing Si + kSt under the constraint that the 
average strain 



E = n{s) = n j dsP{s)a 



(20) 



Pia)=P^{a) 



-l'E(<j) 



(18) 



where the average {■..)ip is taken over the thermo- 
dynamic fluctuations of the dynamic subsystem. 
We thus see that the statistics involves two ingredi- 
ents: the "energetic" , E{(j) and the entropic ST{cr) 
or P,p{<T). On the other hand, all thermodynamic 
quantities (like St = (S't(o'))) are also affected by 
these impacts. 



is known (e.g. measured experimentally). It can 
be estimated as the sample volume dilatation (or 
c-axis expansion for layered compounds). 
According to eq. (18) we arrive at 



P{s) = 



-{•^' — Kct)e 



(21) 



Determining 7' from the constraint (20) we obtain 



P{e) 



-e/E 



(22) 



It should be noted that exactly the same distribu- 
tion would be recovered if only Sj is used as the 
entropy measure. This is a consequence of what 
we called homogeneity, that is the average E is a 
good estimation of the actual strain e, which al- 
most does not fluctuate. Then the experimentally 
available information (20) is relevant to the ther- 
modynamic functions of interest. In other words, 
if the experimental data allow one to recover an 
"exhaustive" picture of the system, then any re- 
finement of the entropy measure (like introducing 
Si + kSt instead of Si) gives no new information. 

Based on the probability distribution the guest 
thermodynamics can be predicted using the exper- 
imental data on the concentration dependence of 
the strain E = E{x) 



St = ST{(i) + aE 



(23) 



Having determined St, we can recover all the ther- 
modynamic functions in the standard way. The 
free energy is also linear in the strain (for simplic- 
ity wc take the same coefficient a): F = F{0)-\-aE. 
Based on this we obtain the guest chemical poten- 
tial as 



dE da 
^ = ^(0) + a— +E— 



dx 



dx 



(24) 



This results coincides with the one obtained in 
our earlier studies [14,15] combining the lattice 
gas model with the continuum elasticity theory. 
The latter implies that the stress can be defined as 
a = dF/dE. If the stress a is composition inde- 
pendent, then from (24) we recover the well-known 
result of Larche and Cahn [18] . 



B. Heterogeneous matrix 

In reality the host matrices are not homoge- 
neous. Even for materials which arc conventionally 
regarded as well-characterized, the heterogeneity is 
practically unavoidable in the experimental condi- 
tions (e.g. because of defects, domain structures, 
etc). On the other hand the amorphous or porous 
matrices are essentially heterogeneous and their 
properties arc known only statistically. In this case 
it is difficult to " invent" a simple experimental test 
(like (20)), giving an exhaustive description. We 
assume that the insertion of guest species into a 
heterogeneous matrix results in a strain distribu- 
tion for various matrix domains. For simplicity wc 
do not consider spatial strain correlations. This 
means that the domain strains Ej are randomly 
distributed according to P{{ei}) = Yii^i^i) 
dependently of the domain position i. 



From the experimental point of view (e.g. mea- 
suring the matrix strain E or from another suit- 
able test) we know only the range in which the 
individual Si can vary. The range of ei varia- 
tion should correlate with the average E. For in- 
stance, if -E = 0, then are symmetrically dis- 
tributed around 0. For concreteness we choose: 
< £i < gi{E). Then, maximizing Si under the 
normalization condition for P{£i), we obtain the 
following step-wise distribution 



Po{e^) 



H{ei)H{gi - e,) 



(25) 



where H{x) is the Heaviside step function. The 
distribution Po{ei) is the one wc can recover using 
the experimentally available information. 

The guest thermodynamics can be predicted as 
before 



ST = 5T(0)+^a,(ei)o 



(26) 



(27) 



where (£,)o = ft/ 2 and ai is assumed to be com- 
position independent. Note that the results above 

are based on the estimation of the internal strains 
(£i)o, using the experimentally available informa- 
tion on the "external" (observable from outside) 
strain _E. As a simple approximation we may sup- 
pose gi{E) = f{Si)E/n, that is, each (ei)o con- 
tributes to the observed strain E proportionally to 
the domain size Si, with J2i{^i)o = E. 

These results can be refined if we explicitly take 
into the host-guest coupling and determine P{si) 
using Si + kSt- This leads to the distribution 



J dsiC'^^i^' 





(28) 



which is more informative in the sense that now 
not all values of are equally probable. The av- 
erage 



l + e""'3'(Ka,ft -1) 
^Q,.(e««i9i - 1) 



(29) 



is different from (£j)o obtained above. It is instruc- 
tive to analyze the origin of this difference. For this 
purpose we expand (si) in terms of the coupling k 



where 



(£») = (£i)o + aiKA° + .. 



A°=[{e^)o-{ei)l]=gf/12 



(30) 



describes the strain fluctuations. The latter are 
related to the domain size fluctuations and thus 



can be called the structural fluctuations. In the 
case of homogeneous matrices all the domains arc 
equivalent, then the structural fluctuations vanish 
A? = and (£j) = (£,)o- Therefore, for heteroge- 
neous matrices our estimation of the strain statis- 
tics becomes more sensitive to structural fluctua- 
tions if we change the entropy measure. This re- 
flects the fact that for heterogeneous matrices the 
actual strain E is different from the one {E) mea- 
sured as the sample dilatation. Their difference 

E-E = ^[(£,) - {e,)o] = J2 (31) 

i i 

takes into account the internal deformations which 
are not detectable experimentally. This result is 
quite general, it does not rely upon the simple 
form for Po{£i) and P{ei) discussed here as exam- 
ples. The only essential fact is that these two are 
different (compare (13) and (18)). The difference 
has its physical origin in the mutual coupling of the 
guest thermodynamics and the host structure. The 
coupling makes it possible to "restore" (at least 
partially) the information inaccessible from exper- 
imental tests. 

Based on the refined distribiition (28) we can 
calculate the guest thermodynamics. 

St = St{0) + J2 + «^ a^AO + ... (32) 

i i 

Considering St as a measure of the guest ther- 
modynamic fluctuations (see above), we conclude 
that the latter are coupled to the host structural 
fluctuations. This reflects the fact that the matrix 
response to the guest insertion is not only due to 
the host nature but also involves the host- guest 
coupling. Comparing the chemical potential 

/x = M(0) + ^ai^(£i)o + «5Za?^A° (33) 

i i 

with (27), obtained using the experimental infor- 
mation, we see that (33) gives a more precise esti- 
mation, but also contains essential physics. In par- 
ticular, if X]i(^«)o Oj then from (27) one would 
conclude that the insertion is almost topotactic 
and that the internal host distortion is irrelevant. 
In contrast, eq. (33) implies that although the ma- 
trix is not strained on average, the expanded and 
contracted domains do contribute to the guest en- 
ergetics. This is observed experimentally [5] on 
various disordered host matrices, which exhibit 
negligible volume expansion in comparison to well- 
structured materials. Nevertheless, the electro- 
chemical response of disordered insertion systems 
is remarkably different. This can be understood 
in terms of two effects. First is the configura- 
tional transitions [14] described by /i(0). Struc- 
tural disorder induces a broad distribution of bind- 
ing energies on different matrix sites. This leads to 



quite steep insertion isotherms with multiple short 

plateaus. A simple model capable of recovering 
this effect can be found in ref. [19]. 

The second effect is due to the structural tran- 
sitions, which arc supplemented by structural fluc- 
tuations (see cq. (33)). In order to analyze this 
contribution within our simple model, let us recall 
that (ei)o = 9i{E{x))/2 and AO = g^{E{x))/12. 
Assuming the simplest non-trivial concentration 
dependence gi{E{x)) = biX, we obtain from (33): 

II = At(0) + 6ii + Kx (34) 

where S/i = J^i^i'^i/'^ induces a shift of the 
isotherm in comparison to the strain free state, and 
K = plays the role of additional re- 

pulsive {K > 0) interaction. Since the chemical 
diffusion coefficient is proportional to dfi/dx, it is 
clear that the structural ffuctuations tend to in- 
crease diffusivity in comparison to the strain free 
state. This seems to be a general tendency for dis- 
ordered matrices, as it follows from experimental 
studies [5,20]. 



V. CONCLUSION 

The coupling of structural and thermodynamic 
fluctuations in the course of various-type insertion 
processes is investigated within a combination of 
Gibbsian statistics and the information theory ap- 
proach. The main focus is on complex (hetero- 
geneous) insertion systems, where the microscopic 
information is not sufficient for making a link be- 
tween the system mechanics (Hamiltonian) and its 
statistics (probability distribution). Then the sys- 
tem can be viewed as an entropic coupling of a dy- 
namic and a stochastic subsystem. It is shown that 
the coupling makes it possible to restore (at least 
partially) the information, inaccessible from exper- 
imental tests. This enables one to make physically 
reasonable predictions under limited information 
on the system. 

It should be emphasized that our results are also 
applicable to matrices which are initially homo- 
geneous but become heterogeneous in the course 
of insertion. For instance, this takes place in Li- 
graphite intercalation compounds, which exhibit 
staging phenomena. In this case the occupied 
and empty graphite galleries are not equivalent. 
Also, if a gallery is partially occupied, then one 
deals with distorted domains [4]. In this case the 
strain fluctuations as well as their spatial corre- 
lation should be taken into account. In a simple 
approximation this physics is reflected by our re- 
sult (see eq. (33)). Our approach is also capable 
of describing the swelling of electroactive polymers 
in the course of intercalation. We plan to address 
these issues in a future study. 
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